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Abstract During April and May 1979 multidisciplinary 
observations were conducted in the central Arctic Ocean from 
three manned ice stations to explore the nature of the submarine | 
Lomonosov Ridge (Weber, 1979). Positions of the LOREX (Lo- | 
monosov Ridge Experiment) ice camps were monitored contin- 
uously using geodetic satellite doppler receivers and the U.S. 
Navy Navigation Satellite System. Near real-time doppler data 
reduction on the ice using broadcast satellite ephemeris provided 
operational horizontal positions with accuracies of about + 250 
m. A precise geodetic satellite doppler reduction program has. 
been adapted to accommodate linear station motion in a simul- 


[ 
: 


*Contribution of the Earth Physics Branch No. 1036; Lorex Contribution No. 14. 


eg 


Marine Geodesy Vol. 7, Number 1-4 
0149-0419/83/030171-00$02.00/0 


BLE. 


BOREAL INSTITULB 
171 ~—~«LIBRARY 


eo Ager 
546770 


‘te eR 


POLARPAM 


ta 
<i 
aa) 
Oo 
a 


o 
pan 9 
<2) 
<a S| 


} 


(LO€x) PL°9ZE" Lggsweg 


> 


dd 


e*. 


+ 


%, 


r 


172 J. Popelar and J. Kouba 


taneous multistation, three-dimensional adjustment in phases 
using postfitted precise satellite ephemeris, if available. Com- 
plete reprocessing of the LOREX doppler data set has produced 
average errors of + 48 m for single pass solution and + 24 m 
for mean three-four horizontal station positions. The station el- 
lipsoidal heights show strongly correlated variations in excess of 
5 m, with the average error of + 0.45 m. Station velocities are 
also strongly correlated, with the mean of about 225 m/h and 
the maximum of 1,240 m/h. Total strain and hourly strain rate 
components of sea ice have been evaluated for a homogeneous 
two-dimensional strain model. The configuration of the ice 
camps provides a unique solution for mesoscale strain, indicat- 
ing a major change in the pattern of pack ice deformation over 
the Lomonosov Ridge. The total strain reflects plastic defor- 
mation, which takes place mainly in episodic events. The strain 
rate also indicates nearly diurnal oscillations, with amplitudes 
an order of magnitude smaller than those corresponding to the 
major deformation events. 


Introduction 


The Lomonosov Ridge Experiment code, named LOREX 79, was 
organized by the Earth Physics Branch in cooperation with the 
Polar Continental Shelf Project. This multidisciplinary expedition 
to explore the nature and origin of the major submarine mountain 
ridge bisecting the Arctic Ocean between the Canadian and Si- 
berian continental shelves took place in the spring of 1979 when 
geophysicists, geologists, and oceanographers carried out mea- 
surements from the drifting sea ice, took samples, and studied the 
ocean floor in the vicinity of the geographical North Pole. Three 
permanent ice camps—the Base Camp, Snowsnake, and Iceman— 
forming approximately an isosceles triangle with sides of about 
60-km long and a 100-km base served as operation bases during 
the experiment (Weber, 1979). The ice camps were established up- 
stream of the narrowest portion of the Lomonosov Ridge, and it 
was expected that the Transpolar Current would carry them across 
the ridge with an average speed of about 5 km per day. With 
known irregularities of sea ice drift it was essential to the success 
of the expedition to provide reliable all-weather positioning of the 
ice camps with the highest achievable accuracy. The Navy Navi- 
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gation Satellite System, also known as Transit, was the only sys- 
tem capable of providing real-time navigation to support field 
operations at the ice camps with a potential for further substantial 
improvement of the positional accuracy using postfitted precise 
satellite ephemeris and highly refined doppler data reduction soft- 
ware. 

Although the Transit system facilitates accurate navigation of 
vessels on high seas and precise geodetic positioning of points on 
the earth’s surface at latitudes up to 70°, the polar orbits of the 
Transit satellites give rise to two undesirable effects in the vicinity 
of the geographic poles: (1) every orbit of every satellite of the 
system is visible, producing at times severe mutual interference 
degrading the data quality; and (2) all the satellites pass nearly 
overhead, resulting in poor fix geometry and, moreover, frequent 
loss of signal at closest approach due to the tracking antennae gain 
pattern. These effects are particularly severe during unattended 
operation of standard Transit satellite receivers. 

Satellite positions of the LOREX Base Camp were evaluated in 
real time with the satellite receiver on line to a disc-based mini- 
computer system most of the time. Raw satellite data from Snows- 
nake and Iceman camps were recorded on magnetic cassettes, 
which were then processed on the main camp computer with a 
three- to seven-day delay. The two-dimensional satellite fix eval- 
uation software developed for BIONAV has been used to process 
the data on the ice during the expedition (Wells, 1976; Wells and 
Grant, 1977). At the same time copies of the raw and majority 
voted data files have been preserved on seven-track magnetic tapes 
for further analysis and postprocessing (Popelar et al., 1981). 

A precise geodetic satellite positioning software package GEO- 
DOP (Kouba and Boal, 1976) has been modified extensively to 
accommodate sequential simultaneous positioning of slowly mov- 
ing stations. The program employs higher order modeling of en- 
vironmenal and instrumental effects, reflects time and space 
correlation of model parameters during simultaneous tracking 
from several stations, and takes advantage of the availability of 
postfitted precise satellite orbits. Thus, much improved positions 
of the LOREX camps have been obtained in three dimensions to 
satisfy the precision requirements of the plumb-line deflection ex- 
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periment (Johnson, 1982). The improved positions have made it 
possible to evaluate differential motion of the three stations and 
provide a unique solution for two-dimensional strain of the sea 
ice during the LOREX expedition. The results are particularly val- 
uable, since for the first time variations of ellipsoidal height for 
the ice stations can be correlated with horizontal sea ice strain 
components. Analysis of the vertical coordinate provides infor- 
mation on variations of sea surface topography in the central Arc- 
tic Ocean. 


Satellite Dynamic Positioning 


Although the real-time satellite navigation satisfied well the op- 
erational requirements on the ice and provided positions accept- 
able for most of the scientific projects carried out during the 
LOREX expedition, the plumb-line deflection experiment re- 
quired satellite positions to about + 30 m to match the expected 
accuracy of astro fixes. High-precision satellite doppler positions 
are routinely obtained by multipass multistation solutions, based 
on postfitted satellite ephemeris, at 1 to 2 m level for geodetic 
purposes. However, geodetic models assume earth-fixed stations, 
whereas motion of an ice station can be appreciable and has to be 
modeled accordingly. 

The new dynamic doppler reduction program GERDOP uses 
comprehensive error modeling of instrumental, orbital, and en- 
vironmental effects for simultaneous multistation three-dimen- 
sional adjustment in phases (Kouba, 1979, 1981). Several program 
modules have been modified or completely rewritten to accom- 
modate station position and velocity variations under the as- 
sumption of linear station motion over specified equal-length time 
intervals. Selection of the time interval is crucial; it has to be short 
enough to permit a linear approximation of the true motion and 
provide for accumulation of sufficient number of satellite passes 
to determine mean station positions and linear velocities. The sta- 
tion motion is modeled by three independent linear fits to succes- 
sive station latitude, longitude, and height coordinates weighted 
according to their formal errors. A new linear fit is initiated at 
the beginning of each interval, and departures from the least- 
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squares line are modeled as exponentially correlated time proc- 
esses in both the station coordinates and velocities to reflect the 
strong positive correlation of station positions over short time in- 
tervals. The doppler data reduction model accepts only dual-fre- 
quency data, corrects for higher order ionospheric effects, and 
respects the time and space correlation of the tropospheric re- 
fraction scaling factors and correlation between doppler data and 
model parameters during simultaneous tracking from several sta- 
tions. Reference orbit bias is obtained from the simultaneous so- 
lution for each satellite pass, with parameter weights providing 
desirable constraints. A sequential adjustment procedure is used 
to update station positions, velocities, model bias parameters, and 
corresponding variance-covariance matrices after each pass. In 
this way only new observations are processed to update results of 
previous solutions without the loss of continuity, thus greatly in- 
creasing the economy of computation. Linear extrapolation is used 
to initialize interval mean station positions and velocities, using 
time- and space-dependent exponential correlation functions. 
Postprocessing of the LOREX satellite doppler data is based on 
majority-voted output obtained during the real-time navigation 
preprocessing. Table 1 gives statistics for the LOREX dynamic 
positioning. The total of 4,920 passes from the three stations that 
have passed the preprocessing data quality checks represents 2,700 
passes of one of the five operational satellites. This means that 
each accepted satellite pass was tracked at an average by two of 
the three LOREX stations, thus effectively taking advantage of 
the simultaneous multistation solution. The very low pass rejec- 
tion rate increases by 20% the overall data recovery in comparison 
with the real-time navigation. Orbital data for satellites 12 and 20 
(1,064 passes) have been derived from the broadcast message in a 
similar way as for the navigation processing, whereas precise post- 
fitted orbits provided by the U.S. Defense Mapping Agency have 
been used for satellites 13, 14, and 19 (1,636 passes). This has 
significant impact on the point-positioning accuracy: the predicted 
broadcast orbits provide satellite positions in + 20-30 m range, 
whereas the accuracy of the precise ephemeris is + 2 m. The se- 
quential adjustment of consecutive passes for which ground tracks 
change direction rapidly in the vicinity of the geographic pole 
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Table 1 
Summary of LOREX dynamic positioning. 

Main Camp Snowsnake Iceman Total 

Station ID SO Sl $2 
Passes tracked cL 2,193 1,940 7,445 
computed 2,189 1,368 1,363 4,920 
accepted 2,162 1,334 1,340 4,836 
Accepted per day Sez. 29.9 30.5 97.6 
Data recovery 64% 62% 69% 65% 
Pass rejection rate 1.2% 2.5% 1.7% 1.7% 


practically eliminates the effect of larger cross-track errors pro- 
duced by the overhead geometry. Moreover, the simultaneous 
multistation solution increases the accuracy of relative station po- 
sitions, which is significant for evaluation of the differential sta- 
tion motion and studies of pack ice dynamics. | 

Mean polar stereographic coordinates of the LOREX stations 
as obtained by the dynamic positioning software for three-hour 
intervals are shown in Figure 1. The overall average error of the 
mean positions is + 24 m in latitude and longitude and + 0.42 
m in ellipsoidal height. These estimates are somewhat conserva- 
tive, since the first ten days of doppler data (day 92 to day 102) 
show much higher noise due to an initial oscillator instability pro- 
ducing much greater formal errors; the height solution is most 
sensitive and had to be practically fixed for that period. It is cor- 
rect to say that errors in mean positions vary from interval to in- 
terval depending on the amount of data, pass geometry, and 
character of ice motion; errors under + 15 m are encountered 
about 40% of the time. 

Figure 2 displays the station velocities and their directions rec- 
tions with respect to the Greenwich meridian. It is apparent that 
motion of the three stations is strongly correlated, reflecting rel- 
ative rigidity of the sea ice cover at scales of 100 km. The average 
station velocity during LOREX was about 225 m/h reaching a 
maximum of 1,240 m/h on day 124. The direction of motion 
changed rapidly at times showing rather large scatter for low ve- 
locities, but the prevailing direction was about 190° with respect 
to the Greenwich meridian. Accumulated track length for the Main 
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Camp was 304.6 km, whereas Snowsnake and Iceman track lengths 
were 241.4 km and 240.6 km, respectively. The direct separations 
between the starting and terminal points for the three stations are 
149.8 km, 140.6 km, and 158.6 km. The LOREX pass summary 
file is available on magnetic tape and lists receiver frequency off- 
set, cartesian and geodetic coordinates, velocity, and separation 
between stations from simultaneous tracking for the 4,836 indi- 
vidual pass station solutions. 

The standard errors for latitude, longitude, height, and station 
separation are derived directly from the solution vari- 
ance-covariance matrix, which also reflects in addition to data 
quality and pass geometry the constraints imposed by the time- 
and space-dependent exponential correlation functions used for 
modeling physically realistic conditions from pass to pass. The 
mean errors for the pass solutions are about + 48 m for latitude 
and longitude, + 0.45 m for height, and + 43 m for station sep- 
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Figure 1. Mean three-hour LOREX ice camp positions. SO = Main Camp, S1 
= Snowsnake, S2 = Iceman. 
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Figure 2. LOREX station velocity plots. 
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aration. These figures are again somewhat conservative due to the 
inclusion of the noisy data from the first ten days of the expedi- 
tion, but they agree well with the errors given above for the three- 
hour mean positions. The slightly smaller error for station sepa- 
ration provides an independent check, as the accuracy of the rel- 
ative station positions is free of orbit-related biases. The station 
velocities and the Greenwich azimuths of the motion have been 
evaluated as a three-pass running average based on consecutive 
station pass solutions. They show a higher mean value of about 
270 m/h compared to 225 m/h derived from the station mean po- 
sitions. This very likely reflects the short-period variation in ice 
motion, which is probably a major factor limiting the achievable 
accuracy of the results. 


Ellipsoidal Heights 


The simultaneous three-dimensional solution provides geometrical 
heights of the electrical centers of receiving antennae above the 
reference ellipsoid. The height time series (Figure 3) show a re- 
markable correlation of station height variations from day 109 on. 
These have amplitudes from 0.2 to about 4 m and range from one- 
day events to systematic trends clearly dominating any changes in 
height differences between the stations. Figure 3 does not contain 
any correction for antenna height, and it seems to reflect clearly 
the fact that the antenna at the Main Camp (SO) was located on 
the roof of the navigation hut about 2.30 m above the ice surface, 
whereas the Iceman antenna (S2) was placed on a small transport 
box only about 0.45 m high. Since it had not been anticipated that 
station ellipsoidal heights could be evaluated from satellite data, 
the antenna heights above water were not measured in the field. 
Approximate heights above ice have been determined from pho- 
tographs of the antennae taken in the field, and it is estimated 
that the uncertainty of the height above water is less than 0.3 m. 

The maximum sea surface slope and its Greenwich azimuth have 
been evaluated from the corrected mean three-hour ellipsoidal 
heights of the LOREX stations (Figure 4). This regional sea-sur- 
face slope with respect to the reference ellipsoid shows rather small 
variations between 5’’ and 15’’, with the mean Greenwich azi- 
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Figure 3. Mean three-hour ellipsoidal heights of LOREX ice camps. 


muth of about 37°. However, such a slope would produce a gen- 
eral reduction of ellipsoidal heights for the LOREX stations along 
their paths, which is inconsistent with the results showing signif- 
icant increase after day 121 (Figure 3). The discrepancy indicates 
the presence of a large time-varying component of the computed 
ellipsoidal heights common to the three stations. There are two 
possible sources of such variations: (1) the first relates to varia- 
tions of the instantaneous sea ice surface resulting from dynamics 
of the ocean, atmosphere, and pack ice; and (2) the second would 
reflect systematic effects due to satellite orbital biases and data- 
reduction model deficiencies. While the former group would be 
of considerable interest to oceanographic studies, the latter rep- 
resents disturbances that cannot be eliminated easily in the present 
data set. However, satellite-related effects are considered small, 
since precise orbits have been used for three of the five satellites, 
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Figure 4. Regional sea-surface slope with respect to the reference ellipsoid. 


and the error-modeling procedures have produced consistently 
high precision results for routine geodetic satellite surveys in sim- 
ilar circumstances. Only the possibility of highly anomalous io- 
nospheric conditions in the polar region may be of some concern; 
it would be therefore highly desirable to include one or two land 
stations in future adjustments of sea ice stations and properly con- 
strain the height solution. 
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Many short-term variations in ellipsoidal height can be corre- 
lated with changes in horizontal station velocities (Figure 2) and 
the horizontal pack ice strain discussed below. The sudden in- 
crease of the station ellipsoidal heights starting on day 121 cor- 
relates with a sharp decrease of the surface water salinity and a 
small increase in water temperature reported at Iceman by Poun- 
der (1980). 


Differential Station Motion and Sea Ice 
Deformation Model 


It is apparent from the shapes of the LOREX triangle as shown 
for day 109 and 145 in Figure 5 that considerable deformation and 
rotation have taken place. The nature of the deformation is re- 
vealed by changes in relative station positions as given by lengths 
and orientation of the sides of the triangle (Figure 6). The side 
connecting the Main Camp and Snowsnake (SO-S1) shows little 
change in station separation before day 118 and after day 127, 
while nearly 2-km extension followed by about 4-km contraction 
took place in between; it is of interest that the side was approxi- 
mately parallel with and crossing the Lomonosov Ridge at the 
time. The side also shows rather steady clockwise rotation of more 
than 15°. The other two sides clearly indicate episodic extension 
taking place in the direction roughly perpendicular to the Lo- 
monosov Ridge with a generally smaller rotation. 

Sea ice deformation includes opening of leads, formation of 
pressure ridges, and shear displacement along cracks and leads. 
These approximately linear features display locally large velocity 
variations and discontinuities. However, on a sufficiently large 
scale the pack ice sheet can be treated as a continuum, assuming 
large numbers of floe-to-floe discontinuities with no single one 
excessively long (Maykut et al., 1972; Coon et al., 1974; Thorn- 
dike and Cheung, 1977). The pack ice sheet deformations are eval- 
uated in two-dimensions due to the great difference between 
horizontal scales (~ 1,000 km) and vertical scales (~ 10 m), and 
because horizontal velocities are basically constant throughout the 
thickness of any particular ice floe. The polar stereographic co- 
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Figure 5. LOREX triangle deformation due to differential sea ice drift between 
days 109 and 145 of 1979. 


ordinates X-northing and Y-easting have been used for the strain 
analysis. 

The following mathematical model has been adopted to express 
displacements v,, and v,; of station j in the X and Y coordinates: 


OV. OV. 
Vxj = ae jXi + ae iJi + X oj 
Ox Oy 
(1) 
dv, dv, 
vice apc ay OS Pisum corte Sp Aim ae b* 
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Figure 6. Relative changes of the LOREX triangle side lengths and orientation. 
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where 
V5 = X(t) — X,(0) 
Xj = X0) — X;(0) 
yi = Y\(0) — Y,(0) 
Le] = 1, 2,..., n number of stations 


[(X,(t), Y,(t)] and [X,(0), Y,(O)] represent the current and 
initial positions of station j 

dv, /0xX, dV,/OxX, dV,/dy, dv,/dy are corresponding displace- 
ment gradients 

X,;, Yj are constants representing the translation 


Equation (1) can be expressed as 


Vs=er,+¢; (2) 


where 


v, is the displacement vector, r; is the position vector, c, is the 
translation component of the displacement vector and e represents 
the deformation tensor matrix. Thus, according to Equation (2), 
e serves as a transformation matrix from the position vector to 
the displacement vector and represents a distortion and/or rota- 
tion of a body. Elements of the deformation tensor and the trans- 
lation components are obtained from the two sets of linear 
Equations (1), providing coordinates and displacements of at least 
three suitably distributed stations are known. If more than three 
stations are involved, a least-squares solution can be obtained as- 
suming homogeneous strain in the region under consideration. 

The nontranslational deformation tensor e defined by the dis- 
placement gradients can be expressed as the sum of a symmetrical 
and an asymmetrical tensor 


e=ert+ao 
with elements 
Ey = 42 (Cx + Cur) = Ex 


Oy = (Cy — Gi) = —On 
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The first part defines the symmetric strain tensor 
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The e,, and e,, strain components represent extension along the X 
and Y axes, whereas 2e,, gives the simple shear. Characteristics 
values (the eigen values) of the strain tensor e define the maximum 
and minimum strain, also called principal strains and their direc- 
tions. These coincide with axes of the strain ellipse representing 
the irrotational part of the nontranslational two-dimensional de- 
formation. Strain invariants are given by 


d ae, naps (5) 


called dilatation (or divergence) and 
En ea [(e,, is Ey)? + 4e2,] 4 (6) 


representing maximum shear. 

The above theory holds for infinitesimal strain characterized by 
small displacement gradients for which the squares and products 
are negligible. However, progressive pack ice deformation shows 
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considerable displacement gradients even for incremental strain 
over time intervals of a few hours. Therefore, the general linear 
transformation 


xe a b| X 
» = ca] é (7) 


is used to relate the final position (x’, y') of a particle after de- 
formation to its initial position (x, y). It is easy to show that the 
strain components can be expressed in terms of the transformation 
matrix elements as follows: 


é =a—l 
€é,, = 2(b+C) 
(8) 
éy=d-1 
w = “’%(b—-—c) 


The finite homogeneous strain components are then evaluated 
from the following expressions: 


{S 
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The mean three-hour station positions have been used for eval- 
uating strain components. The polar stereographic station coor- 
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dinates for day 105 at 15/ UT define the initial ‘‘undeformed’”’ 
triangle configuration. At that moment reliable continuous sat- 
ellite doppler data became available from all three LOREX sta- 
tions. The strain monitoring continued until day 145 at 15/ UT 
when the Snowsnake station satellite tracking terminated. 


Progressive Sea Ice Deformation—Total Strain 


The total strain components for the LOREX triangle have been 
obtained by comparing station configurations with the initial ‘‘un- 
deformed’’ state, using the expressions (7) to (9) for finite ho- 
mogeneous strain. Figure 7 shows the nonrotational strain 
components with respect to the polar stereographic axes. Direc- 
tion of the maximum strain and the principal strains are shown 
in Figure 8, whereas the pure rotation, maximum shear, and dil- 
atation are plotted in Figure 9. The progressive pack ice and de- 
formation is well defined from the LOREX triangle, and two 
distinct deformation patterns can be identified. The first between 
day 110 and 120 is characterized by relatively small deformations 
with direction of the maximum strain nearly parallel with the 
Greenwich meridian. The second pattern from day 127 on displays 
progressive steplike extension in the direction approximately per- 
pendicular to the Lomonosov Ridge with increasing dilatation. The 
two transition periods (before day 110 and between day 120 and 
127) show rapidly changing direction of the axis of the maximum 
strain, and while the first transition indicates significant dilatation 
the second displays just the opposite, i.e., divergence of about the 
same magnitude. Systematic rotation by several degrees in op- 
posite directions is also observed during the two transitions. Day 
121 marks the most dramatic changes, and interestingly enough it 
coincides with the breakup of the Main Camp that was crossing 
the crest of the Lomonosov Ridge at the time. Moreover, as men- 
tioned above, a significant drop in surface water salinity was re- 
corded at Iceman, and a major increase of sea surface heights oc- 
curred at all three stations (Figure 3). These facts can be interpreted 
in terms of distinctly different hydrographic regimes in the Ca- 
nadian and Eurasian basins adjacent to the Lomonosov Ridge 
(Aagaard, 1981). However, one cannot rule out a transient event 
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Figure 7. Nonrotational total strain components. 
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Figure 8. Principal total strain components and their orientation. 
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Figure 9. Pure rotation, maximum shear, and dilatation for total deformation. 
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as a possible source no matter how unlikely the coincidence with 
crossing of the ridge might be. Only additional new data could 
resolve the ambiguity. 

In general, the total deformation during the 40-day period is 
considerable and reaches nearly 30% for maximum strain, exceeds 
8% for dilatation, and shows pure rotation up to 5.5°. 


Incremental Sea Ice Deformation — Strain Rate 


The strain rate calculations are based on differences of consecu- 
tive mean three-hour station positions. Since the magnitude of dif- 
ferential strain effects for three-hour intervals is only marginally 
larger than the uncertainty in the mean station positions and both 
can vary a great deal, it is important to consider the known formal 
position errors when evaluating the strain rate components. This 
has been done by independent least-squares cubic spline smooth- 
ing of the elements of the transformation matrix in Equation (7), 
with weights reflecting the formal errors of station positions. The 
incremental strain components computed from the expressions (8) 
and (9) have been scaled to represent strain rates per hour, and 
these are shown in Figure 10 for the nonrotational components, 
whereas Figures 11 and 12 show the principal strain rates, hourly 
changes in dilatation, maximum shear, and pure rotation. 

The major events reflecting large-scale permanent deformations 
clearly defined by the total strain (Figures 8 and 9) are expressed 
by conspicuous narrow anomalies of the strain-rate components, 
which permit more accurate determination of their timing and 
character. Often rapid, large changes in one direction are followed 
by somewhat lesser change in the opposite direction. Under close 
examination regular, nearly diurnal variations of the strain rate 
components can be observed between the major events. Most no- 
ticeable are these variations in the dilatation rate between day 112 
and 120. Some correlation between strain rate components and 
the fine submeter ellipsoidal height variations (Figure 3) also can 
be observed. Detailed strain and tilt observations on a single large 
ice floe should be considered in the future to determine if this 
oscillation is related to elastic ice deformation due to tidal forces. 
In any case, additional observations are necessary to establish the 
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Figure 10. Nonrotational strain rate per hour (units of 10°). 
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Figure 11. Hourly rates of pure rotation (in degrees), maximum shear, and dil- 
atation (units of 10°). 
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Figure 12. Principal strain rates per hour (units of 10°) and their orientation. 
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suitability of the two-dimensional homogeneous strain model and 
evaluate the effect of large discontinuities between the floes on 
the computed strain components. 

The anomalous hourly changes have reached about 0.3% for 
maximum strain and dilatation, whereas maximum hourly change 
in pure rotation has been 0.15°. The above-mentioned nearly diur- 
nal variations of the same components show amplitudes about an 
order of magnitude smaller. The greatest variance is displayed by 
the direction of maximum strain, which shows signs of discontin- 
uous systematic trends, particularly in the second half of the ex- 
periment. An analysis of statistical frequency produces two peaks: 
one around —65° and one somewhat broader and smaller around 
—10°. The average hourly changes over the duration of the ex- 
periment amount to 0.024% extension for maximum strain and 
0.011% dilatation, which is consistent with the total deformation 
results. 


Conclusions 


The Navy Navigation Satellite System facilitates reliable real-time 
navigation for drifting ice stations in the vicinity of the geographic 
North Pole. The positional accuracy can be improved significantly 
using postfitted precise satellite ephemeris and the newly devel- 
oped dynamic doppler data reduction software. The new software 
accommodates multistation three-dimensional sequential posi- 
tioning of slowly moving platforms under the assumption of linear 
motion over specified equal-length time intervals. It reflects time 
and space correlation of model parameters, including environ- 
mental and instrumental effects, station coordinates, and mean 
interval velocities. The sequential adjustment scheme updates all 
model parameters and corresponding variance-covariance mat- 
rices after each satellite pass, assuring data manageability and ex- 
ceptional economy of computations. In this way all 2,700 satellite 
passes have been reduced in one computer run, using three-hour 
intervals to approximate station drift by linear motion. The over- 
all average error of the mean three-hour station positions is less 
than + 25 m in horizontal coordinates and about + .4 m in el- 
lipsoidal height, while the average station velocity was about 225 
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m/h with the maximum of 1,240 m/h. The average error of the 
4,836 individual pass-station solutions derived directly from the 
variance-covariance matrices is less than + 50 m in horizontal 
coordinates and + .45 m in height, probably reflecting the high 
frequency ice jerk. 

The positional accuracy is sufficient to monitor progressive pack 
ice deformation using a two-dimensional homogeneous strain 
model. The configuration of the three LOREX ice camps provides 
a unique solution for total mesoscale strain, indicating that a ma- 
jor change in the ice deformation pattern takes place above the 
Lomonosov Ridge in the vicinity of the North Pole. It is char- 
acterized by a change in the direction of maximum strain of more 
than 50° and clear divergence when crossing the ridge. The total 
strain reflects mostly plastic deformation, which takes place mainly 
in episodic events when the ice yields under the effects of shear 
and vertical stress. 

Analysis of strain rate can be enhanced by smoothing of ob- 
servations, which takes into account formal errors of individual 
mean station positions. The strain rate components produce con- 
spicuous anomalies corresponding to the major deformation 
events. Moreover, nearly diurnal oscillations with an order of 
magnitude smaller amplitude are quite apparent. However, their 
interpretation would have to include detailed analysis of the suit- 
ability of a two-dimensional homogeneous strain model for pack 
ice deformation studies, which can only be based on presently un- 
available redundant observations. 

Similarly, the interpretation of ellipsoidal heights would require 
at least one land-based station to be included with the ice stations 
in the simultaneous adjustment to constrain the height solution. 
Additional monitoring of oceanographic and meteorological pa- 
rameters is necessary to identify sources of pack ice motion and 
deformation. Future expeditions could provide more information 
by optimizing satellite-observing strategy to assure sufficient re- 
dundancy of measurements and by obtaining complementary geo- 
detic, oceanographic, and meteorological data to study the 
interactions between the ocean and atmosphere, which determine 
ice conditions. This would also be beneficial to other scientific 
projects and contribute greatly to the safety of operations. 
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